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ABSTRACT 



Flow visualization of low speed flow has been obtained 
in this thesis by the injection of a thin stream of gas. 

This injected gas has a refractive index which differs from 
that of the uniform flow in the wind tunnel. A frozen fringe 
hologram of the introduction of this gas indicates regions 
of fringe shift. A pathline is produced by connecting points 
of maximum fringe shift. Pathlines generated in such a man- 
ner around a cylinder agree with the stream functions of 
potential flow theory. 

A rapid means has been developed of data reduction of 
fringe information from a photograph of a reconstructed holo- 
graphic interf erogram. Data points are gathered using a 
digitizer and a mini-computer. The array of data points so 
obtained can be examined and corrected if necessary and fin- 
ally expanded by means of Lagrangian interpolation. The 
expanded array is transferred to a fortran program which ac- 
complishes the integral inversion and subsequently prints 
out the density field. By using this method of data reduc- 
tion and inversion, a large time savings is realized. 
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I. INTRODUCTION 



This thesis consists of two parts. The first part is an 
investigation of low speed flow visualization by means of 
holographic interferometry. The lack of significant density 
changes over a profile at low speeds has made interferometric 
investigations of low speed flow difficult. In this study, 
flow visualization at low speeds has been obtained by inject- 
ing a thin stream of gas in the wind tunnel. This gas had a 
refractive index which varied from that of the uniform flow 
in the wind tunnel. A frozen fringe hologram of the flow 
field with the introduced gas indicated regions of fringe 
shift. The pathline can be found by connecting points of 
maximum fringe shift. Pathlines generated in such a manner 
were compared with potential flow theory by superimposing a 
family of stream functions over the experimental pathline. 

This comparison was done with a cylinder profile in a uniform 
steady flow. 

The second part of this thesis concerns itself with a 
rapid means of data reduction of the large amount of fringe 
data that is obtained in holographic interferometry. A dig- 
itizer and a mini-computer were used to effect this rapid 
means of data reduction. The number of data points is ex- 
panded in the mini-computer by using Lagrangian interpolation. 
The whole array of points are also automatically transferred 
to a fortran data reduction program. This existing fortran 
program then inverts the integral equations to determine the 
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density field. In the past, fringe data had been collected 
from a hologram by the use of an optical comparator to deter- 
mine the amount of fringe shift along each plane. The data 
was typed onto IBM cards and fed into a fortran program for 
reduction. The present development of the use of a digitizer 
and the ability to automatically feed data into a fortran 
program greatly reduces the time to obtain a density field 
from a hologram. 
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II. FLOW VISUALIZATION 



A. HOLOGRAPHY AND FLOW VISUALIZATION 

Dennis Gabor in 1948 (1) was the first to invent a meth- 
od of wavefront reconstruction called holography. With the 
development in the early 1960 ' s of a coherent source of ra- 
diation, the laser, the field of holography rapidly developed. 

Holography is the process of stoi'ing wavefronts in a 
photosensitive material called a hologram. At a later time 
these wavefronts may be reconstructed from the hologram to 
form an image which has the same depth and parallax charac- 
teristics as the real object. This is accomplished by re- 
cording the interference pattern between the light diffracted 
by the object and a reference beam of light. An optical 
arrangement is set up to split the coherent light source into 
a reference beam and a scene or object beam. Within certain 
restrictions the paths of these two beams are arbitrary. 

The photosensitive plate is placed in such a position that 
the reference beam and the diffracted light from the object 
fall upon it. As stated in Ref. 1, each point on the holo- 
gram contains the full information of the object as seen 
from that point. Full information constitutes intensity and 
phase information. The phase information is derived from 
the interference pattern. When the object is reconstructed, 
each of these points superimposes its own complete picture, 
which has depth and parallax. The recording of phase infor- 
mation (with the associated dc h and parallax) i. what 
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distinguishes holography from photography which only records 
intensities . 

A frozen fringe hologram is one which has been doubly 
exposed. One exposure records the reference situation. The 
second exposure records whatever perturbation has been intro- 
duced. In the present low speed flow visualization experi- 
ment the reference condition of exposure is the steady flow 
of uniform gas over the two-dimensional profile. After the 
reference exposure an optical component is moved slightly to 
provide fringe lines. The flow field is perturbed by the 
introduction of a thin stream of gas into the uniform steady 
flow ahead of the profile. A second exposure is taken on 
the same hologram plate of the perturbed flow. In the re- 
construction of the resulting hologram, the two recorded 
flow fields interfere with each other and show a deviation 
in the superimposed finite fringes. Since both exposures 
utilize the same optical components, the fringe lines are 
lines of constant differences in optical path lengths. Any 
imperfection in the optical quality of the apparatus is thus 
eliminated . 

When a gas of a different density (optical path length) 
is carefully introduced into the flow field, it is swept 
through the wind tunnel and follows a pathline. Holography 
is extremely sensitive to variations in optical path length 
and by using holographic techniques a line of fringe shifts 
can be seen in the region of the injected gas. By carefully 
adjusting the amount of gas introduced in relation to the 
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wind tunnel velocity, a small well defined line of fringe 
shifts can be produced. A pathline can be created by con- 
necting points of maximum fringe shift. 

B. EXPERIMENTAL APPARATUS 

The arrangement of the experiment is shown in figure 1(a). 
This type of set-up is known as a modified Z-Type Schlieren. 

A Korad Q-switched laser (ruby wave length 6943 A) is mounted 
on an optical rail. A small five-milliwatt helium-neon con- 
tinuous gas laser is placed adjacent to the rail to provide 
alignment of the optics. The beam is expanded and projected 
onto a parabolic mirror. This mirror collimates the beam to 
facilitate easier reconstruction. The reference beam travels 
through the beam splitter and is reflected off two mirrors 
before illuminating the holographic plate. The object beam 
is reflected off the beam splitter and then projected through 
the wind tunnel. Just before the tunnel the object beam 
passes through a diffusing plate and a grid. The actual ex- 
posure time with the Q-switched laser was on the order of 
twenty nano-seconds. This stopped all movement during the 
exposure, but did not prevent movement between exposures. 

The tunnel was mounted separately from all other components 
to minimize unwanted movement caused by vibration. Alignment 
of the optics was done by projecting the small laser beam 
through the ruby laser head and out into the optical system. 
The initial alignment of the small laser with the Q-switched 
laser was accomplished with burn holes in photographic paper. 
This technique is explained in Appendix A of Ref. 2. Fringe 
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lines were formed by moving the beam splitter .003 inches be- 
tween exposures. The beam splitter was mounted on a preci- 
sion X-Y translating stage. 

The wind tunnel is 3.5 inches by 3.5 inches and construct- 
ed out of plexiglass (not of optical quality). The speed of 
the tunnel can be varied from zero to twenty-two feet per 
second. A section of honeycomb was placed at the entrance to 
serve as a flow straightener and as a support for the tube 
which injects the gas. 

The tube used for injecting the gas had an outside dia- 
meter of 3/32 inches. It was connected by tubing to a flow 
meter. This meter allowed visual adjustment of the injection 
gas flow rate. 

Two profiles were used, both constructed from plexiglass. 
Pictures of the profiles and their supports are shown in Fig. 
1(e). A cylinder was used having a length of three inches 
and a diameter of one inch. The second profile was a wing 
having a length of three inches, a chord of 1.5 inches and a 
thickness of 3/16 inches. Both profiles were supported by 
two upright pieces of plexiglass. A grid was lightly inked 
onto the supports to provide a reference during reconstruc- 
tion. The ends of the profiles were painted black for better 
definition of their outline. There was a snug fit between 
the pins, which extended from each end of the profiles, into 
the two supports. This fit allowed the wing to be set at a 
desired angle of attack and remain there. When the profile 
was adjusted between the supports then the assembly was slid 



into the tunnel. The honeycomb was placed at the entrance 
and the injection tube inserted through the honeycomb. 

Reconstruction was done with a fifteen-milliwatt helium- 
neon laser. The beam was expanded from .04 inches to .25 
inches by the use of two mirrors. A sketch of the recon- 
struction is shown in Fig. 1(d). A shutter was placed be- 
tween the second mirror and the hologram plate holder. In 
order to photograph the real image the camera lens was re- 
moved. A dark cloth was used to cover the laser and prevent 
any extraneous light from overexposing the film. 

Other equipment used in this experiment: 

1. power supply and circulator for Q-switched laser 

2. eye protectors 

3. Agfa 10E75 holographic emulsion plate and Polar- 
oid Type 55 film 

4. auto collimator for internal Q-switched laser 
alignment . 

C. EXECUTION AND RECONSTRUCTION 

Before taking a frozen fringe hologram, the alignment of 
the optical system was verified using the small 5-milliwatt 
laser. While this was being done the power supply for the 
Q-switched laser was turned on. The recommended warm-up 
time for the power supply is five minutes. The circulator 
was also turned on to stabilize the laser head and the eta] on 
at a low temperature. A small amount of air was blown onto 
the optical surfaces to remove dust. The mirror directing 
the alignment beam through the ruby was removed and a 2 . 5min 
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transverse mode selector pin hole was placed in its holder. 
The Agfa 10E75 emulsion plate was placed in the plate holder 
which was under a dark cloth. After appropriate warnings to 
other personnel in the area and donning safety glasses, the 
laser was charged. The wind tunnel was set at the desired 
speed and allowed to stabilize. The black cloth covering 
the holder was removed, the laser fired, and the black cloth 
replaced. Then the beam splitter was moved .003 inches. 

The injection gas was turned on and the flow regulated. The 
capacitor bank for the laser was again charged and, after 
removing the cloth, fired. The laser, wind tunnel and the 
injection gas were turned off while developing the emulsion 
plate . 

Emulsion plate development was standard: 

1. Kodak D-19 developer for five minutes 

2. acetic acid stop bath for thirty seconds 

3. fixer for five minutes. 

The plate was washed in water followed by a rinse in Kodak 
Photoflo wetting agent. A hot air dryer was used to dry the 
plate. Occasionally the plate was overexposed and a bleach 
was used. Whenever this was done an index matching solution 
of xylene was also used to reduce the noise level of the 
hologram. Further information on bleaching and index match- 
ing is available in section 10. S. 5 of Ref. 1. 

Reconstruction was done with a fifteen-milliwatt Stabi- 
lite laser using the setup explained in section B and 
illustrated in Fig. 1(d). The virtual imago was quickly 
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reconstructed to assess the quality of the hologram. This 
was easily done by using the collimating lens of the Stabi- 
lite laser and looking through the hologram. If the quality 
of the plate warranted a photographic record then the plate 
was reversed in its holder. The collimating lens was re- 
moved from the laser and the beam was slightly expanded to 
reduce speckle. The virtual image was projected onto the 
frosted glass at the back of a lensless camera. For a dis- 
cussion of real and virtual images the reader is directed to 
Ref. 1. An f-16 aperature was used on the shuttering device. 
Black cloth covered the' camera, holder, shutter and the la- 
ser. The plate holder was adjusted to present the correct 
picture on the ground glass of the camera. A photographic 
record was made using Polaroid Type 55 positive/negative 
film. 

D. COMPARISON OF RESULTS WITH POTENTIAL FLOW THEORY 

Two-dimensional potential flow theory was used to give a 
comparison of experimental pathlines obtained by connecting 
points of maximum fringe shift. A discussion of this theory 
is presented in Ref. 4. 

The program listed in Appendix A calculates stream func- 
tions of uniform flow past a circular cylinder with circula- 
tion. The stream function ip is given as: 

* = Uy — + f- In (— ) 

r J r^ 2u a 

where U = uniform velocity ; r =Vx 2 + y 2 

a = cylinder radius; t = circulation. 



15 



An example of the graphical output of a cylinder with circu- 
lation is shown in Fig. 2(c). In the experimentation for 
this thesis there was no circulation around the cylinder. 

When the photograph of the cylinder in Fig. 2(a) is super- 
imposed over the stream functions calculated from the program 
in Appendix A then Fig. 2(d) is the result. Notice the close 
agreement between the stream function labeled 3 and the 
points of maximum fringe shift over the profile. Downstream 
of the cylinder the shift in the stream function 3 is be- 
lieved due to the heavier carbon dioxide used as the injected 

gas. As the wind tunnel speed is increased the point of di- 
vergence of the injected gas with the stream function is 
shifted further downstream. It was found that for very low 
tunnel speeds the fall of the carbon dioxide occurred at 
earlier points. At wind tunnel speeds of 22 feet per second 
the line of fringe shifts caused by injecting carbon dioxide 
agreed favorably with a stream function over a cylinder. 

E. DISCUSSION OF FLOW VISUALIZATION 

Three types of injection gases were investigated. Helium, 
argon, and carbon dioxide were introduced into the flow field. 

Helium and carbon dioxide provided an acceptable fringe shift. 

The most advantageous gas was carbon dioxide since its high 
molecular weight prevented the injected stream from rapidly 
diffusing in the tunnel. Fig. 2(a) shows that this fringe 
shift is small and well defined over the profile. 
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To accurately define the pathline the injected gas must 
be swept through the tunnel. If the flow rate out of the 
injection tube is too high, then the subsequent path of the 
gas does not follow a pathline. To ensure pathline accuracy 
the velocity of injection gas is matched as closely as pos- 
sible to the velocity of the tunnel. Both the tunnel veloc- 
ity and the injected gas velocity were 22 feet per second. 
Velocity matching required calibration of both the wind tun- 
nel and the injection flow metering device. Once the veloc- 
ities are matched then the amount of gas injected is dependent 
on the injection tube diameter. At low speeds a tube which 
had a diameter 5% of the height of the tunnel was able to 
define a fringe shift throughout the hologram. At higher 
tunnel speeds a smaller tube can be used since the gas flow 
rate is increased and more gas is introduced to cause a 
shift. Fig. 2(b) is an example where the flow velocities 
were not matched. Notice in the figure how the injected gas 
is blown onto the lower surface of the wing. The Kutta Con- 
dition is easily recognized in the circled area at the 
trailing edge. 

F. CONCLUSION DERIVED FROM FLOW VISUALIZATION 

The sensitivity of holography to slight variations in 
density makes this flow visualization technique ideal for 
low speed flow investigation. Depending on the nature in 
which you inject the gas one can emphasize different features 
of the flow. Although it was not done for this thesis, 
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perhaps vortices can be visualized by injecting more gas. 

As shown by Fig. 2(a), when the velocity of the injected gas 



is matched to that 
fringe shift makes 



of the wind tunnel the resultant line of 
a pathline visible around the profile. 
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III. DATA REDUCTION 



A. RAPID DATA REDUCTION TECHNIQUE 

Due to the large amount of data generated by holographic 
interferometry, a rapid means of data reduction is desired 
in the determination of a density field from holography. To 
date the emphasis has been on developing a better fortran 
program which inverts the integral equations in a calculation 
of a density field. A comparison of the different methods 
used to invert the integral has been discussed by Sweeney 
[Ref. 6]. This paper addresses the problem of data acquisi- 
tion and transfer of these data to a fortran inversion pro- 
gram. A digitizer is used in determining the fringe shift 
in an interf erogram obtained by reconstructing the frozen 
fringe hologram. A mini-computer processes the data using 
Lagrangian interpolation, and then stores it. At an appro- 
priate time the data are recalled from the storage tape cas- 
sette in the mini-computer and automatically transferred to 
a large computer for density calculation. In the previous 
method of data reduction the fringe data had been collected 
by the use of an optical comparator which determined the 
amount of fringe shift along each plane of a hologram. The 
fringe data were then typed onto IBM cards and submitted to 
a fortran program in a large computer. The present develop- 
ment of the use of a digitizer and the ability to store data 
and then automatically transfer .it to a fortran program 
greatly reduces the time needed to obtain a density field. 



To obtain a photograph of fringe data, a doubly exposed 
hologram is taken of the aerodynamic phenomenon. The first 
exposure is of the reference condition which usually is a 
no-flow situation. Flow is then introduced into the area in 
question and a second holographic exposure is taken. Finite 
fringe lines are established on the hologram by a slight 
movement between exposures of an optical component. The 
holographic plate is developed and then the scene is recon- 
structed using a laser source. A photographic record can be 
taken of the real or the virtual image. For a further ex- 
planation of holographic interferometry the reader is direct- 
ed to Ref. 5. 

The actual integral inversion in the fortran program can 
be done in many ways. A comparison of different methods has 
been discussed by Sweeney [Ref. 6], Some of these methods 
are : 

1. Representation of the integral by orthogonal poly- 
nomials 

2. Fourier Expansions 

3. Fourier Transforms 

4. Fast Fourier Transforms. 

The overwhelming amount of data generated by holographic 
interferometry requires a rapid method of data acquisition , 
evaluation, storage, and transfer. 
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B. EQUIPMENT DESCRIPTION AND OPERATION 

This data reduction technique uses a Hewlett Packard dig- 
itizer (HP9864A), calculator (HP9830A) and printer (HP9866A). 
A photograph of the apparatus is shown in Fig. 3. An acous- 
tic coupler was also used to allow the calculator to act as 
an on-line terminal with the Naval Postgraduate School IBM 
360 computer. The power supply for the digitizer was mounted 
under the table and the printer, as can be seen in Fig. 3, 
is resting on top of the calculator. 

Data were introduced into this equipment with the use of 
the cursor attached to the digitizer. A button is on the 
cursor which defines the origin on the magnetic table. The 
coordinates of single points are entered by depressing but- 
ton "S". Continuous points can be entered by depressing "C" . 
This last feature was not used in the program. 

The programmable calculator uses Basic language and has 
an initial memory capability of 1760 words. Additional mem- 
ory can be added externally or internally to the calculator. 

A tape cassette with 32,000 word storage capacity can be 
used with the calculator. Program or data can be stored on 
or retrieved from this cassette. The entire program is on 
the first portion of this cassette. The remainder of the 
magnetic tape is used to store information for up to nine 
pictures, each picture containing twenty fringe lines and 
twenty planes. By employing a Terminal Read-Only-Memory 
(Term ROM) the calculator becomes an on-line terminal for an 
IBM 360 computer. While in the Term ROM mode the calculator 
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retains its basic operating ability. It is possible, al- 
though not advisable, for an operator to execute a terminal 
command and while the IBM operation is being performed, use 
the calculator locally. Operation of the calculator as a 
terminal was in accordance with the normal terminal procedure. 

C . PROCEDURE 

The procedure of using this method of data reduction is 
extremely straightforward. A listing of the programs to 
transfer data is given in Appendix B. Each of the six pro- 
grams is stored in a file on a tape cassette. Each program 
is referred to by the cassette file in which the program is 
stored. The arrangement of the programs were designed to 
allow an operator to collect and expand fringe data for up 
to nine pictures. At a later time, when the terminal mode 
for the IBM 360 is available, then the stored data are sent 
to the large computer for reduction. A flow diagram of 
these operations is shown in Fig. 4. A description of the 
procedure used in the data reduction follows. 

A photograph of a frozen fringe hologram was taped to the 
table of the digitizer. If the hologram did not have a well 
defined grid than a thin transparent grid was taped over the 
photograph to provide a reference. The calculator, printer 
and digitizer were turned on. The data reduction tape cas- 
sette was placed in the calculator and the command "LOAD" 
was executed. The program was started at point A on Fig. 

4 by executing the command "RUN". Statement 20 requested 
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the operator to enter the time, day, month and year. These 
numbers were introduced into the calculator by typing them 
onto the display. All numbers were separated with commas. 

The display was executed and the program (HP file 0) contin- 
ued to run. HP file 0 printed the time, day, month, year, 
and the program organization. At statement 190, HP file 1 
was automatically retrieved from the cassette and loaded in- 
to the calculator. HP file 1 asked for the orientation of 
the grid. The origin and an arbitrary second point on the 
horizontal axis were entered by depressing the cursor "0" 

(for the origin) and then the cursor "S". By lifting the 
cursor i inch off of the table the established origin was 
erased and the orientation procedure had to be repeated. 
Whenever the calculator is ready to accept a data point from 
the digitizer a soft beep is heard. The calculator asks for 
the number of fringe lines to be used from the picture. Next 
the number of planes in the picture is entered. The program 
is able to accept up to 20 fringe lines and 20 planes. The 
equation of each fringe line is entered by entering two arbi- 
trary data points on the line in the linear region. After 
all lines were entered the program prints the mean and stand- 
ard deviation of the slopes of all lines. Also the differ- 
ences from the mean of each slope are printed. The operator 
makes an evaluation as to the validity of his fringe lines 
by comparing the slopes. If any slope varied greatly from 
the mean then all fringe lines can be reentered. An example 
of the printout is shown in Appendix C. If the fringe lines 
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are satisfactory then "CONT" was executed and HP file 2 is 
loaded into the calculator. The picture number and the num- 
ber of data points in each plane are entered. As each data 
point is entered the elevation and fringe shift are printed. 
At conclusion of each plane the array of elevations and 
shifts are stored on a cassette tape file. The first number 
in the array is the number of data points in that plane. 

The second number in the array is the elevation of the first 
data point. The third number is the fringe shift for the 
first data point. Elevations and shift alternate through 
the remainder of the array. 

The reference file is filled after acquiring all the data 
points in the picture. The reference file must be compatible 
with the fortran data reduction program. A program based on 
Fourier Transform methods, [Ref. 3], was modified to accept 
data input parameters in a specific order. After the refer- 
ence file is filled, it is stored in the cassette under the 
appropriate picture number. The parameters in the reference 
file are shown in Appendix B, file 2 starting with statement 
430. The operator was at point B in Fig. 4. Either he 
can gather more data from different pictures or expand the 
data for fortran reduction. 

HP file 4 is loaded for data expansion. Expansion is 
done using Lagrangian interpolation. All of the data points 
from the plane to be reduced are printed out. The program 
asks for the data point numbers which will be the basis of 
the array of expanded data. The number of points to use for 
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curve fitting should be less than seven. If a greater number 
of points are used then the computational time greatly in- 
creases but more importantly the curve fit to the data points 
is erratic. If more than one picture needs to be reduced 
then the same point numbers are automatically used for data 
expansion of each picture. The first number in the expander 
array is the number of data points to follow. The second 
number is the number of the plane in the picture. Before the 
expansion of the array all positions were initialized to 
9.99999. This is shown in File FT03F001 of Appendix E where 
the actual data that has been fed in are indicated and those 
positions in the array that are not expanded data are 9.99999. 
The reason for initializing the array is for convenience of 
data transmission. The absence of the nines following the 
array should be investigated by reading out the stored data. 
After the expanded data points are stored in the cassette 
the operator loads HP file 3. 

File 3 is the program which automatically feeds the data 
into the fortran system. Directions for its use are printed 
by executing "LIST 10,110". To reduce the transmission time 
the expanded data is sent to the IBM 360 in arrays of six 
numbers. The arrays of six numbers are placed on a storage 
disk in the IBM 360 system under a file labeled "FILE FT03F001" . 
The fortran reduction program calls FILE FT03F001 for its 
data input. 
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D. DISCUSSION 



The previously described technique permits extremely 
rapid data acquisition and transfer. The ability to rapidly 
acquire, transfer, and store fringe data in large quantities 
allows the operator to process many holograms. To determine 
a density field around an asymmetric object in a wind tunnel, 
as many holographic views as possible are used. Each of 
these holograms is taken at a different angle and provides 
fringe data as seen from that angle. This program has the 
ability to process nine holograms of 400 data points each. 

An improvement of this data reduction method can be made 
by incorporating a plotter peripheral and by developing a 
better interpolating routine. At present the expanded data 
is based on six fringe shift data points. When more than 
six points are used the expanded array becomes erratic. One 
would like to utilize as many of the available 20 data points 
as possible on each plane. A new interpolating routine could 
be developed to accept more than six points and still yield 
accurate results. The expanded array and the original data 
points could be compared on the plotter. From the graphical 
output of the plotter an improved subset of data points could 
be chosen to expand the array. By incorporating this scheme 
not only the best order polynomial is used to expand the data 
but also the best set of data points. 

With the ability to evaluate and store up to nine sets of 
fringe data on a cassette, valuable .operator/computer time is 
not lost transferring and correcting fortran input data. 
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E. CONCLUSION 



A procedure for rapid fringe data collection, transfer 
and storage has been accomplished. With slight modification, 
this procedure is applicable to all fortran data reduction 
techniques. With one photograph of 20 fringe lines and 20 
points, the use of a digitizer was 20 times faster than that 
of an optical comparator. This was assuming the optical com- 
parator was used in conjunction with a programmed desk calcu- 
lator to calculate fringe lines and fringe shifts. The amount 
of time in the transfer of data to the fortran program was 
15 times less by using this technique. When many pictures 
are taken to completely determine the density field around 
an object, the magnitude of time saved is very large. 

The ability to rapidly store and transfer the large 
amounts of data obtained from a finite frozen fringe holo- 
gram results in a considerable time savings. 
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Figure 1(a). Schematic of holographic arrangement. 
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Figure 1(b). An overall view of the holographic experiment 




Figure 1(c). Cylindrical profile mounted between 
plexiglass supports 




Figure 1(d). 



Profile and supports 
Notice honeycomb and 



inserted in tunnel, 
injection tube. 
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Laser 




Figure 1(e). Schematic of reconstruction 
of a hologram. 
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Figure 2(a). Photograph of a reconstructed hologram 
of a cylinder in a uniform flow of 
22 fps. Notice the line of fringe 
shifts over the profile caused by the 
injection of a thin stream of carbon 
dioxide . 
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Figure 2(b). Photograph of a reconstructed hologram 
of a two-dimensional wing in a uniform 
flow of 6 fps. Notice the Kutta Condi- 
tion at the trailing edge of the profile. 
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Figure 2(c). Graphical output of potential flow program 
given in Appendix A. The stream functions 
are for a uniform velocity of 10 ft /sec 
around a cylinder of 1.18 in. with a clock- 
wise circulation of 75 ft 2 /soc. 
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Figure 2(d). Fringe lines of Fig. 2(a) superimposed on 

graphical output of stream functions around 
a cylinder. 




Figure 3. A view of the Hewlett Packard digitizer, 

printer and calculator. On the digitizing 
table is a finite frozen fringe hologram 
with a transparent grid and the electronic 
cursor. To the rear of the table is the 
acoustic coupler. 
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I data to be expanded and reduced, see Fig. 4(b) 



Figure 4(a). Flow diagram of data transfer program. Dashed 
lines are manual operations. Program file num- 
ber and statement number is in the right upper 
corner of each operation. 



37 




TAPE III: Expanded fringe 

shift stored for each 
picture 



recall reference file 
and — — — — 
expanded data 



I 



n 

* 



file 3 , *■ 
Transfer Data 
for Data 
Reduction 



-10 



IBM 360 
Reduce Fringe 
Shift Data to 
Obtain a Density 
Field 



-> 



Density 

Field 

Printed 




Figure 4(b). Flow diagram of the data transfer program. 

During the file 3 operation the calculator 
is acting as a terminal . 
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APPENDIX A 



Listing of Potential Flow Program 
This program calculates stream functions around a cylinder 
with circulation. 



10 


REM*** 


* 2 — IT 


FLOW** 


12 


DIM fiC 


250 J 




20 


DISP " 


ENTER 


XMINj 


30 


INPUT 


X 1 ? X 2 


j Y 1 ? Y2 


40 


SC RLE 


X 1 J 72 


? Y 1 ? Y 2 


50 


XflX I S 


0 5 + 1 




68 


YflXIS 


0 J + 1 





70 REM FOR TITLE TYPE "LETTER" 

8 0 R E M P 0 S I T 1 0 N !•! 3' T H 13 1 0 Pi .R Y R E Y S 
90 REM "EXECUTE" IS RETURN/LINE FEED 
1 0 0 R E M " S T O P " T E R M 3. N fl T E S I... E T T E R M 0 13 E 
110 LFi E3 E L < * j 4 ? 2 ? 0 ? 2 / 3 > 

120 DISP " TYPE TITLE" 

130 LETTER 

i 40 D I SP " I NCR Y < , 3 ) , Z < < 250 > » PH I CO NV < . 05 > " 5 
150 INPUT Y6 ? Z ? Z 1 

1.60 DISP "ENTER PICTURE DIfi&RCT VELOCITY" 5 

170 INPUT R'ZjUI 

175 Y=0 . 5 

180 U=R2*U1 

190 R=R2/2 

200 T=-4*PI*R*U 

2 1 0 D I S P " C R I T I C fl L. C 1 R C U L fl T 1 0 N I S " T 
228 WRIT 4006 

2 3 8 D I S P " E N T E R C I R C U l. fl T 1 0 N < C W 1 S H EG) " j 
240 INPUT T 



2 5 0 D I S P " F 0 R 0 1 R C L E G 0 C 0 N T 5 9 0 

260 STOP 



70 


PLOT X 1 ? 0 r + 1 




i t nr 
. I* --J 


DISP "POSITIVE 


PHI BELOW X AXIS" 


76 


WRIT 4000 




: s 0 


DISP "ENTER PHI 


j PHI END j PH I I NCR 


:90 


INPUT Pj P2 ? P3 




: 0 0 


PLOT Xlj0j+i 




:08 


1 = 1 




: 0 9 


1 1 -• ( X2--X 3 )/Z 




00 


FOR X=X1 TO X2 


STEP 1.1 


! 2 0 


Y5=Y6 




130 


IF X=0 THEN 350 




:40 


GOSUB 440 




:50 


RC I J=Y 




151 


DISP I 





39 



• j cr ( tr 


1 = 1 + 1 






360 


NEXT X 






361 


1 = 1 






362 


MSP " PH I = " P 






363 


WRIT 500 






365 


FOR X=X1 TO 


v ■ j 


STEP n 


366 


PLOT XfflCI]> 


+ 2 




36? 


1 = 1 + 1 






363 


NEXT X 






370 


PEN 






380 


IF P=P2 THEN 


4 


20 


396 


P--P+P3 






410 


GOTO 300 






420 


STOP 






440 


N=0 






450 


B - 3 L-! R ’• X * X ■ V 


Y > 




468 


fi=IJ*Y - ( IJ*R*R 


# Y 


) •••'’ \ X * X + Y * Y • } * L 0 G b 


465 


fl=-fi 






470 


Q=R-P 






480 


IF RBS(fl-P) 


\ := 


2 J THEN 538 


490 


IF Q*N >= 0 


THEN 510 


500 


Y5=Y5/2 






5 1 0 


IF FKP THEN 


55 


0 


520 


Y=Y+Y5 






5 3 0 


N=Q 






540 


GOTO 450 






550 


Y---Y-Y5 






568 


H=Q 






570 


GOTO 450 






580 


RETURN 






590 


P L. 0 1 -■ R i 0 :> + 1 






600 


FOR X=-R TO 


R 


STEP 0.01 


6 1 0 


Y=SQR<R*R-X* 


X > 




620 


PLOT X? Y j +2 






630 


NEXT X 






t* 4 0 


P L 0 1 R ? 0 j + 2 






650 


FOR X=R TO - 


R 


STEP --0.01 


6 6 0 


Y=-SQR<R*R-X 


-e- x 


j 


670 


PLOT X j Y ? + 2 






680 


NEXT X 






690 


P L U f ~ R j 0 ? + 3 






700 


PEN 






710 


GOTO 270 






720 


END 
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APPENDIX B 



Listing of Calculator Data Transfer Programs 

B.l Listing of HP File 0. This program prints the organiza- 
tion of the files in the tape cassette. By executing "Run 
200" the stored data of any picture may be printed. 



10 REM THIS PROGRAM BY DAV.ID POPE USES fi CALCULATED PRINTER 

1 1 R E n R N D II I G ITIZER C H P 9 8 3 0 R ? H P 9 8 6 6 R ? H P 9 8 6 4 fl ) T 0 R E C 0 R D T H E 

12 REM FRINGE SHIFT FROM Fl PHOTO OF R RECONSTRUCTED HOLOGRAM 

13 REM THE INFORMATION IS FED INTO THE CP /CMS SYSTEM OF RN 

14 REM IBM 368 COMPUTER WHICH DOES THE RCTUflL DATA REDUCTION 

15 REM THE NORMAL. FLOW OF DATA IS: FILE 0?FILE 1 ? FILE 2i 

16 REM FILE 4? AND THEN FILE 3 FOR THE CP /CMS HOOKUP. 

19 FIXED 0 

20 DISP "ENTER TIME DRY MONTH YEAR"* 

30 INPUT Y1 ? Y2? Y3? Y4 

4 8 P R I N T " T I M E " Y 1 " D A Y " Y 2 " M 0 N T H " Y 3 ” Y E A R " Y 4 

5 8 P R I N T " D I R E C T 1 0 N S A N D 0 R G A N I Z A T 1 0 N " 

7 0 prj; nt "ijp TO 29 PLANES WITH 20 POINTS PER PLANE" 

88 PRINT "DATA FOR 9 PICTURES MAY BE GATHERED AND STORED" 

90 PRINT "FILE 0: ORGANIZATION? READOUT 

100 PRINT "FILE 1: EQUATIONS FOR UP TO 28 LINES" 

110 PRINT "FILE 2: FRINGE DISPLACEMENT FROM THE LINES" 

120 PRINT "FILE 3: MATE WITH CP/CMS SYSTEM" 

1 3 8 P R I N T "FILE 4 : LAG R A N G I A N I N T E R P 0 L. A T 1 0 N 

140 PRINT "FILE 5: COMPARES TWO PLANES? REDEFINES DATA FILES" 
168 PRINT "FILE 6: STREAM FUNCTIONS FOR 2-D POTENTIAL FLOW" 
181 PRINT "FILES 7 ON ARE DATA STORAGE FILES" ' 

183 PRINT " STORED DATA READOUT PRESS LOR DO EX?RUN200 EX" 

190 LOAD 1 

2 0 8 R L M * / * * * R E: A D U LJ 1 U L S I 0 R E: ] J D fl ! fl * $ # * 

2 2 0 D I M D S E 4 1 3 ? S I E 5 3 ? T S E 6 3 ? Q S E 1 5 0 3 
238 DISP "ENTER PICTURE NUMBER".? 

235 FIXED O 
246 INPUT P 

2 4 1 P R I N T " P I C T IJ R E N U M B E R " F : ' 

242 FIXED 5 

2 5 U Z 2 — P * 2 3 ~ 1 6 

268 LOAD DATA Z2?S 

2 7 0 L. 0 A r.i D A T A < 2 2 + 1 ) ? T 

2 7 1 L U A D D fl T A '.. Z 2 + 2 ? Q 

273 PRINT "INTERPOLATED DATA FOR CP /CMS" 
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274 PRINT "THIS IS FROM PLANE"QC23 

275 PRINT "POINT #" <Q[ 1 3+3) "SHOULD BE 9. 99808 " 

283 FOR I 1 TO (QC 13+3) 

284 PRINT 1'jQCI] 

285 NEXT I 
296 Z=SC 1 3 

3 8 1 P R I N T " N U M B E R 0 F P L. A N E S ~ " S I 1 3 

302 PRINT "NODE OF OPERATION =="SE23 

303 PRINT "NUMBER OF DATA. POINTS FOR CP/CMS ~" 

304 PRINT "NUMBER OF VIEWS = " SC 4 3 

305 PRINT "SYMMETERY OF FIELD ="SL53 

306 PRINT "REFERENCE DENSITY (MG/CC) =*’TC 1 D 

307 PRINT "WAVELENGTH IN MICRONS ""TE23 

308 PRINT "RMAX (CM) ="TE33 

309 PRINT "FILTER = "TC43 

310 PRINT ”TH (VIEWING ANGLE IN DEG) ="T[53 

311 PRINT "DATA SYMMETERY ="TE63 
378 FOR Zl=l TO Z 

380 Z2=Z1+(P*23)-14 
390 LOAD DATA Z2? D 

4 0 0 P R I N T D [ 1 3 " D A T A P 0 1 N T S I N P L A N E " Z 1 
410 FIXED 3 

420 I "DC 1 3*2 

43U FuR K=2 TO I STEP 2 

440 J=K+ 1 

4 5 0 PRINT " E L E V A T 1 0 N = " D C l< 3 » " S H I F T = "DC ,J 3 
460 NEXT K 
470 NEXT Z1 

7** 1 F' fc I [■-] T ' T*i‘ ■£* "A” *A‘ “jr iv •**;• ■*!;■ ~i" ■+■ "0 “li* “i-- -lj- -*!■■ ■¥■ •+■ 

472 PRINT 

473 PRINT 
480 STOP 
500 END 
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B.2 Listing of HP File 1. This program calculates the equa- 
tions for up to 20 fringe lines in each picture. The mean 
of the slope and the differences from the mean of each slope 
is also printed. Two data points are entered through the 
digitizer to calculate each fringe line. 



1 0 C IJ M Z 3. ? J .1. !> J 2 ? ii b L 2 U ] ? B b L 2 0 .1 j U b L 2 0 J 

19 PRINT "FOR CORRECTIONS YOU MUST START OVER 

20 U ISP "ORIENT THE CURVES" 

30 MBIT 2000 

40 MSP "ENTER THE BEGIN? PRESS. . CONT. EX" 

50 STOP 

60 MSP "ENTER SECOND POINT ON HORIZ AXIS" 

70 NR I T 2O00 
S O N R i 1 E 9 > * ) 

90 ENTER < 9 s * > X j Y 
300 Ji=SIN<RTN<Y/'N>> 

110 J2= : C0S(RTN< Y/N) ) 



120 

130 

140 

150 

160 

170 

180 

190 




30 
40 
260 
270 
“80 
8 1 



DEG 

PRINT "SLOPE OF THE PHOTO IS" <RTN< 
MSP "# OF FRINGE LINES ( II RX OF 28 
INPUT N 

MSP "NUMBER OF P ERNES? MAX OF 29“ 5 
INPUT Z 

PRINT "NUMBER OF LINES "N 
PRINT "NUMBER OF PLRNES"Z 
MM MSI 20] 

E : -0 

FOR Ml TO N 

MSP “FIRST POINT FOR LINE “I 
NR IT 1880 
WRITE < 9 ? * > 

ENTER ( 9 j * > X ? Y 

y, i — x j 2 + Y j i 

Y1=Y*J2-X*J1 

D I S P " S E C 0 N I) P 0 1 N T L I N E " I 
WHIT 1000 



DEGREES 
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cn ; -.n cn ai 



290. WRITE (95*) 

306 ENTER <9»*)X>Y 
3 2 0 X 2 X * J 2 + Y * J 1 
330 Y2=Y*J2-X*J1 
3 4 0 MCI] = ( Y 2 -- Y 1 ) •••- < X 2 - X 1 > 

350 E-E+MC I 3 
360 R=SQR<MC I 342+1 ) 

370 fit I 3=MC I :i.-"R 
380 BC I ]=(--l )/R 
390 CC I 3~'(Y1"T1[ I 3*X1 )/R 
400 Ii ISP "LINE "I" 

410 WAIT 1000 
420 NEXT I 

430 H 1 = E N 

431 FIXED 4 

432 PRINT "THE MEAN OF THE SLOPES IS 



IS ENTERED' 



N 1 



440 S=0 

450 FOR 1=1 TO N 

4 6 0 P R I N T "LINE" I T'DI F F E R E N C E F R 0 1'l M E A N I S 
470 S=S+(Mt I 3-N1 >+2 
480 NEXT I 
490 Ii=SQR(S.--N> 

00 PRINT '"STANDARD DEVIATION OF SLOPES "D 
1 0 D I S P " S A T I S F I E D W I T' H S L 0 P E S ? » . . C 0 N T ,, E X 
30 STOP 



40 LOAD 2 



MC I 3 HI 
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B.3 Listing of HP File 2. For each data point entered this 
program calculates the fringe shift from the appropriate 
fringe line. The array of fringe shifts from each plane are 
stored on tape. A reference file is also filled and stored 
on tape. 



5 C 0 i ' 1 Z I ? J 1 ? J 2 ? fl S C 2 0 ] ? B S [ 2 9 ] j C S C 2 0 3 

is run GSC413 

15 PRINT "CORRECT EACH PLANE BEFORE NOV I NO TO HE NT PLANE" 
20 PRINT "IF POINT CORRECTION THEN, . . STOP CO NT 250" 

2 5 D I S P " C 0 R R E C T 1 0 N . . . S T 0 P C 0 N T 2 5 0 " 

30 WAIT 5800 

35 MSP "ENTER PICTURE NUMBER" 5 
40 INPUT P 

4 5 P R I NT " P I C T U R E N U M B E R ” P 
50 Zl=l 

60 Z2=Z1+ CP+23) - 1 4 

65 FIND Z2 

66 I) ISP "START ENTERING POINTS" 

67 WAIT 1500 

7 0 P R I N T " D I S P L A C E M E N T I N P L A N E " Z 1 

71 PRINT "POINT # ELEVATION SHIFT 

75 MSP "NUMBER OF POINTS IN PLANE " Z 1 ? 

88 INPUT N 



85 GC 1 3=N 



90 J=l. 

95 1 = 1 

100 WRITE <9>*> 

105 ENTER ( 9 ? ) X > V 
125 X1=X*J2+Y*J1 
1 3 0 V 1 Y * J 2 ~ X * J 1 
135 J=J + 1 



140 GC J 3=Y 1 
145 J=J+1 

1 5 0 g c j :i = a c i :i * x 1 + b c i 3 * y i + c c i j 

155 K= J- 1 
160 FIXED 4 

165 REM GC1) IS THE NUMBER OF POINTS 
170 REM G(..J) IS THE SHIFT? G CIO IS THE 
1 7 5 P R I N T I ? T A B 1 0 ? " Y = "GC K II ? T A B 3 6 ? G C J 3 
180 IF I =N THEN 191 
185 1=1+1 



IN THE PLANE 
ELEVATION 



198 GOTO 100 
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191 WRIT 4000 
195 STORE DATA Z2»G 
200 IF Z1=Z THEN 215 
205 Z1=Z1+1 
210 EG TO 60 

215 DISP "SATISFIED WITH THESE POINTS'?" 

220 WRIT 3000 

225 H ISP " CO NT... EX" 

230 STOP 

235 REM EACH FILE CONTAINS ONE PLANE 
240 GOTO 395 

2 4 5 R E M * * * * * C 0 R R E C T 1 0 N S * * * * * 

250 DISP “ENTER LINE NUMBER" 5 
255 INPUT II 
2 6 0 W R I T E ( 9 ? * ) 

265 ENTER <9j*)XjY 
285 X1=X*_J2+Y*J1 
290 Y1=Y*J2-X* Jf 
295 R=2* 1 1 
300 GC R 3=Y 1 
305 R=R+1 

3 1 0 GC R 3= AC 1 1 3*X 1 +BC 1 1 3* Y 1 +CC 1 1 3 
3 1 5 P R 1 N T " C 0 R R E C T 1 0 N P L fi N E " Z 1" L I N E "II? 
3 20 II I S P "IF SR T I S F I E D G 0 . . . C 0 N T 1 9 5 " 
325 WRIT 5080 



330 GOTO 
REM : 



395 



250 

*** FILL 



THE REFERENCE FILE 



400 DIM S I C 5 ] ? T SC 6 3 
415 Z2=P*23-1S 
420 FIND Z2 



•i *-icr 
4.^1 -J 



SC 1 3=2 



430 


DISP 


"ENTER 


MODE OF OPERATION <3} 


435 


I NF'IJT 


SC 2 3 




440 


DISP 


"ENTER 


NUMBER OF VIEWS'! 1 V S 


445 


INPUT 


SC 4 3 




450 


DISP 


"ENTER 


SYMMETRY OF FIELD 1 ! 0} 


455 


I NPUT 


SC 5 3 




465 


D I SP 


"ENTER 


# OF POINTS FOR CF'/'C 


470 


INPUT 


SC 3 3 




505 


STORE 


DATA 


Z z ? y 


5 1 0 


D I SP 


"ENTER 


REFERENCE DENSITY (MG 


5 1 5 


INPUT 


TC 1 3 




520 


DISP 


"ENTER 


WAVELENGTH IN MICRON 


525 


INPUT 


TC 2 3 




530 


D I SP 


"ENTER 


F'iMfiK ( UM " 5 


535 


I NPUT 


TC 3 3 




540 


DISP 


"ENTER 


FILTER"? 


ET £T 

•J fc + -J 


INPUT 


TC 4 3 




550 


DISP 


"ENTER 


THC VIEWING ANGLE IN 


tr tr cr 


INPUT 


TC 5 3 




560 


DISP 


"ENTER 


DATA SYMMETRY (0> " ? 


565 


I NF'IJT 


TC 6 3 





5 






i :r. 



EG 



SHIFT = " GC I 
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,570 


Z2=22+ 1 












575 


STORE 


DATA Z2j T 












580 


PR I NT 


"IF YOU NR NT 


TO 


REDUCE 


THIS 


PICTURE THEN, 


. . LORD 5 


590 


PRINT 


"IF YOU NR NT 


TO 


GATHER 


MORE 


DfiTfl FOR NEXT 


PICTURE 


600 


PR I NT 


ii n 












610 


ni-sp ' 


•DRTfl STORED 


FOR 


PICTURE 


"P 






620 


STOP 














630 


END 
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B.4 Listing of HP File 3. This program recalls the array of 
expanded data from the tape cassette and then the calculator, 
acting as a terminal, automatically sends the data to FILE 
FT03F001 in the IBM computer. 



5 P R I N T " F 0 R D I R E C 1 1 0 H S P R E S S L I S T 1 9 5 1 1 6 . . . E X “ 

19 REM SET ACOUSTIC COUPLER DUPLEX TO 1/2 ft MB TURN ON 
28 REM DIAL 2611 ON PHONE AND PLACE IN COUPLER 
21 REM YOU SHOULD HERR A HIGH PITCHED TONE 

30 REM RED LIGHT ON? THEN THE COUPLER IS ON 

3 1 R E M N H I T E L I G H T 0 N N H E N T R A N S M I T T I N G / R E C I E V I N G 
48 REM TYPE IN "TERM. .. EX" 

50 REM "CP-67 ONLINE" RPPERRS IF CP •••-CMS IS READY 
60 REM TYPE LOGIN PROCEDURE 
78 REM TYPE "I PL CMS" 

88 REM TYPE "TTY FILE FT03F801 78" 

81 REM LAST # IN FT03F . . SHOULD BE PICTURE # 

82 REM SHOULD D ISP" READY PAPER TAPE IN RFADFR. , , " 

83 REM TYPE & TRANSMIT TWO VIRTUAL CARDS NORTH OF TITLE 

90 REM THIS IS THE READ IN DATA FILE FOR THE CP/CMS 

91 REM TYPE IN TITLE <#i> AND TRANSMIT 

92 REM TYPE IN TITLE (#2) AND TRANSMIT 

108 REM PRESS "RUN. . . EX" ? THIS STARTS THE CALCU LATER 
1 1 8 REM .TRANSM I T ANYTH I NG ".CONTROL ( F6 > -S-TRAHSN I T < F9 > 

1 2 0 D I M A I [ 5 I » B S IS]? S S I 1 5 8 I 
125 FIXED 5 

138 DISP "ENTER PICTURE NUMBER " 5 

140 INPUT P 

150 Z2=P*23-16 

160 LOAD DATA 22 5 A 

1 7 0 L 0 A D D A T A < 2 2 + 1 > i B 
i 8 0 L U A D DAT A < Z 2 + 2 ) S 
183 RE 3 I-INT (SC 1 ]> 

1 8 4 P R I N T " Y 0 U A RE RE D U C I N G " A E 3 .1 " D A T A P 0 1 N "f S I N C P / C M S " 
185 PRINT 

190. DISP "ARE YOU READY TO SEND DATA?" 

288 WAIT 3008 



205 


AC 3 3= I NT (SC i :J> 






2 1 0 


D I S P "IF S 0 > P R E S S C 0 N T . . . E X 


i n 




220 


STOP 






236 


N R I T E ( 4 j 3 4 0 > RE 1 3 :■ AC 2 3 > fl C 3 3 ? 


RC 4 I ? R 


C 5 3s 19 


235 


PR I NT AC .13 5 RC 2 ] 5 AC 3 ] • ML 4 ] j AC 


5 3 




240 


NR IT 3000 






258 


NR IT 3000 






260 


W R I T E ( 4 j 3 5 0 > EC 1 ] j BC 2 I , B C 3 1 , 


BC 4 I j B 


E 5 3s BE S 3s 19 


265 


P R I N T B C 1 ] 5 B C 2 ] ■ B C 3 3 5 B C 4 3 5 B C 


5 3s BC 6 


3 


2 6 6 


PRINT 






270 


NR IT 3000 






280 


NR IT 3000 






296 


G=INT<(RC33/6>+l ) 






295 


6=6+1 






300 


FOR 1=1 TO (G*6> STEP 6 






3 i 0 


N R I T E ( 4 > 3 6 0 > 8 E I 3s SC I + 1 3 ? S C I 


i .-i - i !“■■■ r 
"i 1 j ? •::> l. 


I +3 3 > SC 1+4 3 


3 i 5 


P R I N T S Cl]? S C I + 1 ] 5 S C I + 2 3 5 S C I 


+3 J ? SC 


1+4 3? SC 1+5 3 


326 


NR IT 3000 






330 


NR IT 3000 






340 


r I.J R M H T 5 F o .. 0 j B 






350 


FORMAT 6F10. 5? B 






360 


FORMAT 6F 11.5? B 






380 


NEXT I 






390 


N R 1 T fr. 4 ? 3 9 5 ) 1 9 ? 1 9 ? 1 9 ? 1 9 ? 3. 9 j 


i 9 




395 


FORMAT B ? B 5 B ? B ? B ? B 






400 


PRINT “DATA SENT TO FILE IN 


CP. •-CMS 


" 


40 i 


P R I N T " T 0 R E D U C E D fi T R < E X E C U T 


E POPE 


1 > $ POPI 


402 


PR INF " UU 1 PU i 1 N f- I Lb F 1 Bor O 


0 1 " 




4 0 3 


PRINT 






4 1 0 


STOP 






420 


END 







4 0 



B.5 Listing of HP File 4. This program expands the fringe 
shift data by using Lagrangian interpolation. The array of 
expanded fringe data is stored on the tape cassette. 



1 0 D I M X S C 20 3 s V S [ 2 0 ] , D S C 4 1 ] > G S E 1 5 0 3 > U I E 2 0 3 



15 


PRINT 


"TO S 


TOP PRINTING . , 


■ „ . DEL 520 j 520. 


1 6 


PRINT 








17 


WRIT 5 


000 






2 y 


D ISP " 


ENTER 


PICTURE S, PLRi 


-iE • 


30 


INPUT 


Pi ? 22 






40 


33 ISP " 


ENTER 


# OF PICTURES' 


;i 


50 


I NPUT 


Q 






60 


PRINT 


"# OF 


PICTURES TO BE 


■■ REDUCED IS" Q 


70 


H=U=0 








80 


P=R=U[ 


1 3=1 






90 


N=S7=0 








100 


FOR I 


= 1 TO 


150 




110 


GC I 3“ 


9 . 9 9 






120 


NEXT 


i 






1 90 


PRINT 


" P I c 


T U R E " P 1 j " PLANE 1 


: .“j 


200 


—.1 

iLs k!._ u... 


+P1*2 


_ •] A 




2 1 0 


LORD 


DhTR 


?3n T\ 




220 


K=DE 1 3 






230 


J ~ -*' 1 








240 


N~0 ' 








250 


X=DE J 


2 3 






260 


Y=DC J 


2+ i J 






270 


J2=J2 


4 2 






280 


N=FN2 


r 






290 


IF N< 


K THEN 250 





300 IF H-it0 THEN 415 
320 G 03 LIB 660 

330 Ii ISP "# OF POINTS TO USE (<?) 

340 INPUT N 

350 FOR 1=1 TO M 

360 33 ISP "USE POINT #"5 



370 


INPUT 


UC I 3 




374 


FIXED 


0 




375 


PRINT 


"POINT" UE I 3 "USED" 




376 


FIXED 


4 




380 


NEXT I 




390 

400 


33 ISP ' 


'MIN ELEV.jMRX ELEV. 


OF POINTS"? 


4 1 0 


INPUT 


R7 j B7j S8 




4 1 5 


PR I NT 


"INTERPOLATION BASE! 


13 ON "M" DATA POINT 


420 


PR I NT 


"MIN ELEVATION "A7 




430 


PRINT 


"MAX EL.EVATI0N"B7 




440 


PI' I NT 


"ii OF POINTS "S3 





50 



460 6C 2 ]=Z2 
470 11=2 

4 8 0 S 7 = ( B 7 - fl 7 ) • ••" S 8 - 1 ) 

500 FOR X=A7 TO B7 STEP S7 
510 GC 1 1 + 1 ]=FHAX 
515 UISP a 1~1) 

520 PRINT ai-njXjGt 1 1 + 1 ] 

525 GC 13= 1 1-1 
530 11=11+1 
540 NEXT X 

5 y 0 S i U RE D H 1 A P 1 * 2 3— 1 4 ) j G 

591 PRINT "*****“ GC 1 I" POINTS GENERATED AND STORED*#*** 

592 PRINT 
b O O H = H + 1 

610 IF H=Q THEN 640 
620 F'1=P1 + 1 
630 GOTO 80 

6 4 0 D ISP : 'E N D 0 F L A . I N T E R . " 

641 WRIT 4080 

642 PRINT "TO REDUCE DATA LORD 3" 

643 PRINT 
650 STOP 

660 PRINT "POINT #">" X‘V V" 

665 FIXED 4 

670 FOR 1=1 TO M 

680 PRINT I,XC I 3? YE I ] 

690 NEXT I 
695 PRINT 
780 RETURN 
710 END 

720 DEF FN2 < Z ) 

730 FOR I=N TO 1 STEP -1 
740 IF X< XC I 3 THEN 790 
750 NCI + i3=X 
760 YCI+1]=Y 
770 RETURN N+l 
780 XC 1 + 1 ]=XC I ] 

790 YC 1 + 1 3= YE I J 
800 NEXT I 
810 1=0 
820 GtiTO 750 
830 END 

840 DEF FNFKZ) 

8 6 0 b = 0 

980 'FOR 12=1 TO M 
990 T=1 
1000 I=UC 123 
1G10 FOR J 1 = 1 TO M 
1O20 J = UC J 1 3 
1030 GOSUB 1G80 



r.i 



1040 NEXT J1 
1050 S=S+T*YCI] 

1060 NEXT 12 

1070 GOTO 1120 

1080 IF XCI3=XCJ] THEN 1110 

1090 T = T * < Z - X C J ] > / ( X [I]- X C J ] 

1110 RETURN 

1115 T =0 

1120 RETURN S 

1130 END 



B.6 Listing of HP File 5. The beginning of this program 
compares two pictures at one plane. The second section of 
this program redefines all of the data files. 



10 RE M C 0 M P fl R E 2 P L. fl N E S . . . R U N j E X 
2 0 R E M R E D E F I N E D R T fl F I L E 8 . . . R U N 3 4 0 j E X 
30 REM POTENTIAL FLOW OVER 'CYLINDER. . „ RUM570 ? EX 
40 DISP "COMPARE THE TWO PLANES" 

50 DIM ESC 41 3? FSC 41 ] 

60 WAIT 3000 

70 DISP "ENTER THE PLANE NUMBER" 5 
80 INPUT Z 

90 DISP "ENTER PICTURE NUMBERS” ? 

100 INPUT PI 5 P 3 

1 1 0 P R I N T " P L. A N E Z ” P I C T U R E S " P 1 , P 2 

120 Z 1=2+ (Pi *23) -14 

130 Z2=2+<P2*23)-14 

140 LOAD DATA Z 1 ? D 

150 LOAD DATA Z2jE 

160 FIXED 4 

170 IF DC 1 ]<EC 1 3 THEN 230 

180 R=EC 13*2+1 

190 FOR I -2 TO R 

200 FC I 3=DC I I-EC 13 

210 NEXT I 

220 GOTO 270 

230 R=BC 13*2+1 

240 FOR 1=2 TO R 

250 FC I 3=BC I 3--EC I 3 

260 NEXT I 

270 FOR 1=2 TO R STEP 2 
280 J=I+1 

2 9 0 P PINT " B I F F ” F C I 3 , T A B 2 0 > F C J 3 
300 NEXT I 
3 1 0 S 3 0 F‘ 

320 GOTO 48 
339 END 



340 REM ****REDEFINES H FIT ft FILE 
350 FIND 7 
360 FIND ? 

370 FOR P=1 TO 10 
388 22=8*23-16 
390 FIND 22 

400 DISP "PICTURE IS "P 
410 lift IT 500 
420 MARK 1.-5 
430 FIND (22+1 > 

440 DISP "PICTURE "P" ft" 

450 WRIT 200 
460 MARK 1-12 
470 FIND (22+2) 

480 DISP " PICTURE "P" B" 
490 WAIT 200 
500 MARK 1 ? 300 
510 FIND (22+3) 

520 DISP " PICTURE "P" C" 

530 MARK 20-82 
540 NEXT P 
550 STOP 
568 END 
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APPENDIX C 



Printed Results of the Calculator Programs 
C.l Results from File 0, File 1, and File 2. 



TIME 1110 DRY 14 MONTH 5 YERR i974 
li I R E C T 1 0 N S R N D 0 R G R N I 2 fl T 1 0 N 
UP TO 20 PLANES WITH 20 POINTS PER PLANE 
DRTR FOR 9 PICTURES MAY BE GATHERED AND STORED 
F I L E 0 ‘ 0 K G A N I Z A T 1 0 N ? R E A D 0 LJ T ? C 0 M P A R t 2 P I C T U R E S 
FILE 1: EQUATIONS FOR UP TO 26 LINES 
FILE 2: FRINGE DISPLACEMENT FROM THE LINES 
FILE 3: MATE WITH CP/CMS SYSTEM 

F I L. E 4 = L A G R A N b IAN I N T E R P 0 L A T 1 0 N ( N 0 R M A L L. Y N 0 T U S E 3 
F I L E 5: LAG R A N G I A N I N T E R P 0 L. A T 1 0 N U S E D 
FILE 6! COMPARE TWO PICTURES AT A PLANE 
FILES 7 ON ARE DATA STORAGE FILES 
STORED DATA READOUT PRESS LOADS EX ? RUNS 00 EX 
FOR CORRECTIONS YOU MUST START OVER 
S L 0 P E 0 F THE P H 0 T 0 I S 9 . 7 S S 3 9 1 9 44 DEGREE '• 

NUMBER OF LINES 5 
NUMBER OF PLANES 3 



THE MEAN OF 


THE SLOPES 


IS 


-0. 83 


15 




LINE 1.S000 


D IFF ERE NO 


E FROM 


MEAN 


IS 


0 . 004 l 


LINE 2.0000 


DIFFERENC 


E FROM 


MEAN 


IS 


-0 , 0095 


LINE 3.0000 


D IFF ERE NO 


E FROM 


MEAN 


IS 


0 . 0273 


LINE 4.0000 


DIFFERENC 


E FROM 


MEAN 


IS 


-0 . 0330 


LINE 5.0000 


DIFFERENC 


E FROM 


MEAN 


IS 


0.0111 


STANDARD DEV 


I AT I OH OF 


SLOPES 


. Q „ 82 


83 




CORRECT EACH 


PLANE BEFORE MO 


VING 


T 0 


NEXT PLANE 


I F P 0 1 N T C 0 R R E C T 1 0 N T H E N . . . S 


TOP C 


ONT 


2 5 0 


PICTURE NUMBER 3 










DISPLACEMENT 


IN PLANE 


1 








POINT # 


ELEVATION 






SHIFT 


1 . 0800 


Y= 0.31 


5 6 






0. 0085 


2. 0000 


Y- 8.35 


ib? 






0. 0016 


3. 0000 


Y= 0.3S 


! b 3 






0. 0378 


4 . 0000 


Y= 8.41 


r •-» 






0 . 0330 


5. 0800 


Y= 0.43 


• r* id 






0. 0722 


DISPLACEMENT 


IN PLANE 


2 b 0 0 0 0 








POINT 4 


ELEVATION 






SHIFT 


1 . 00011 


Y= 0.25 


i'57 






8. 0537 


2. 0000 


Y= 0.2S 


1 C| O 






8 . 8528 


3. 0000 


Y= 0.2992 






0. 1233 


4 . 0000 


Y = 0.3107 






6. 1 146 



DISPLACEMENT 


IN PLANE 3.0000 




POINT # 


ELEVATION 


SHIFT 


1 . 6000 


Y= 0.2827 


0.01 98 


2. 0000 


Y= 0.3238 


0. 0022 


3. 0000 


Y= 0 . 3650 


0.0559 


4 . 0000 


Y= 0.4044 


--0.0047 


5. 0000 


Y- 0.4258 


0.0677 



IF YOU WANT TO REDUCE THIS PICTURE THEN. .. LOAD' 5 

IF YOU WANT TO GATHER MORE DATA FOR NEXT PICTURE. ..LOAD 1 

T 0 S T 0 P P R I N T I N G .... D E L 5 2 0 » 5 2 0 . „ E X 



# OF PIC 


TURES TO BE REDUCED IS 


i 


PICTURE 


3 PLANE 2 




POINT # 


l iB l 


1 *: 
T 


1 . 0000 


8.2597 


O . 0537 


2 . 0000 


0. 2893 


0 . 0528 


3 . 0000 


8.2992 


O. 1233 


4 . 0000 


0.3107 


0. 1146 


POINT 1 


USED 




POINT 2 


USED 




POINT 3 


USED 




POINT 4 


USED 





INTERPOLATION 


BASED ON 


4.8008 DATA POINTS 


MIN ELEVATION 


0 . 2500 




MAX ELEVATION 


0 . 3200 




# OF POINTS 13 


1 . 0000 




1 . 0808 


0,2500 


0 . 3274 


2. 0000 


8,2558 


0. 1389 


3.0000 . 


0 . 2 b 1 7 


0 . 0223 


4 . 0000 


8.2675 


-0 . 0354 


5 . 0000 


0.2733 


-0 , 0478 


6. 0000 


0.2792 


— 0 , 0254 


7 . 0000 


0, 2850 


0 . 0 1 64 


8 . 0000 


0. 2908 


0. 0657 


9 . 0000 


0.2967 


0. 1094 


10. 0000 


0 . 3025 


0.1 3 4 8 


1 1 . 0000 


0 . 3083 


0. 1289 


12. 8000 


0 . 3 i 42 


0 . 8789 


13. 0000 


0. 3200 


-0 . 9282 


***** 13.0000 


PO I NTS 


G E N E R A T E D A N D S T 0 R E D * * 



TO REDUCE DATA LOAD 3 
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C.2 Printed output of readout of stored data. Listing of 
stored data obtained by executing "Run 200" on File 0. 



PICTURE NUMBER 3 
I N T E R P 0 L ft T E D D A T R F 0 R C P / C M 3 
T H I S I S F R 0 M P LAME 2 . 0 0 0 6 0 
POINT # 16.00000 SHOULD BE 9'i 



1 „ y 0 y 0 0 0 



1 . 0 0 0 0 0 j. 3 . 0 0 0 0 0 

2.00000 2.00000 

3. 00000 0. 32742 

4. 00000 0. 13893 

5 . 0 0 0 0 0 0 . 0 2 r 2 3 0 

6 . 0 0 0 8 0 ~ 0 . 0 3 5 3 9 

7 . 0 8 0 0 0 - 0 .847 0 8 

8 . 0 0 0 0 0 - 0 .025 4 4 

9 . 0 0 0 0 0 0 , 0 1 6 4 2 

1 0 . 0 0 0 0 0 0 . 0 6 5 6 7 

1 1 . 8 8 8 0 0 8 . 1 0 9 4 3 

3 . 2 . 8 0 0 0 0 0 . 1 3 4 8 1 

1 3 „ 8 0 0 0 0 8 . 1 2 8 9 3 

1 4 . 0 8 0 0 0 0 . 0 7 8 8 8 

15. 08800 -0. 02821 

1 6 . 0 8 0 0 0 9 9 9 9 9 9 . 0 8 0 0 0 
HUMBER OF PLANES = 3.08088 

M 0 1) E 0 F 0 P E R A T 1 0 H = 3 . 0 0 8 0 0 

HUMBER OF DATA POINTS FOR CP/CMS = 13.08088 

HUMBER OF VIEWS = 1.08800 

SYMMETERY OF FIELD = 0.68088 

R E F E R E H C E D E N S I T Y ( M G / C C > * 1 . 7 4 5 8 8 

WAVELENGTH IN MICRONS * 0.69448 

RMAK (CM) -- 2.78008 

F I L T E R = 8 . 0 0 0 0 0 

TH (VIEWING ANGLE IN DEG) = 0.88800 
D A T A S Y M M E T E R Y = 8.0 0 6 0 6 

5 . 0 8 8 0 0 D A T A P 0 1 N T S I N P L A N E 1 . 8 8 8 8 0 



ELEVATION-- 0.316 


SHIFT* 


0 . 0 0 8 


ELEVATION = 6.357 


SHIFT* 


0. 002 


ELEVATION- 0.386 


SHIFT* 


0 . 038 


ELEVATION* 8.418 


SHIFT* 


6.033 


ELEVATION* 8.437 


SHIFT* 


8.072 


4.008 DATA POINlS IN 


PLANE 2.008 




ELEVATION* 8.268 


SHIFT* 


8.654 


ELEVATION* 0.289 


SHIFT* 


8.G53 


ELEVATION* 0.299 


SHIFT* 


0. 123 


ELEVATION* 0.311 


SHIFT* 


6. 115 


5.088 DATA POINTS IN 


PLANE 3.008 




ELEVATION* 6.283 


SHIFT* 


0 . 82.0 


ELEVATION* 6.324 


SHIFT* 


6,. 062 


ELEVATION* 6.365 


SHIFT* 


0 . O 5 6 


ELEVATION* 0.464 


SHIFT*- 


6. 685 


ELEVATION* 6.426 


SHIFT* 


O. 068 



■v ■¥- •*!- -r- • 



■V V I -h -'A -V- -P Y -I; Y ■)- j; -l- I, l» ■>*■ -I; “P -I; ■ 
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APPENDIX D 



Fortran Subroutine Which Inputs Data 

This subroutine is in the Fortran data reduction program. 
The title, reference numbers, and fringe shift data are read 
from FILE FT03F001. 



SUBROUTINE DATA ( I T, QA, QB, QS ) 

C ****** THIS INPUTS DATA FROM FILE FT03F??? ****** 

REAL IT(12),QB(6),QS(150) 

INTEGER QA(5) 

READ (3,981)( IT( I ), 1=1,12) 

WR I TE ( 6, 1 00 0 ) ( IT( I ), 1=1,12) 

READ (3,982)(QA( I ), 1=1,5) 

WR I T E ( 6 , 1 0 0 1 ) ( QA ( I ) , 1=1,5) 

READ (3,983)(QB( I ), I =1,6) 

WRITE (6, 1002) (QB( I ), I =1,6) 

I D=6* ( QA( 3 ) / 6+1 ) 

DO 980 I =1. I D. 6 
READ(3,984) (QS( l+J-1), J=l,6) 

WR I TE ( 6 , 1 0 0 3 ) ( QS ( I + J - 1 ) , J = 1 , 6 ) 

1000 FORMAT (IX, 6 A 8) 

1001 FORMAT ( IX, 5 I 8 ) 

10D2 F 0 RM AT ( 1 X , 6 F 1 0 . 5 ) 

1003 FORMAT (1X,6F11.5) 

980 CONTINUE 

981 F ORMAT ( 6A8 ) 

982 FORMAT (518) 

983 FORMAT ( 6F 10 . 5 ) 

9 84 FORMAT ( 6 F 1 1 . 5 ) 

I X=QA ( 3 ) 

DO 985 I =1, I X 
QS ( I ) =QS ( 1+2) 

985 CONTINUE 
RETURN 
END 
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